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Abstract. We present an extension of the methods of classical Lie group analysis 
of differential equations to equations involving generalized functions (in particu- 
lar: distributions). A suitable framework for such a generalization is provided 
by Colombeau's theory of algebras of generalized functions. We show that under 
some mild conditions on the differential equations, symmetries of classical solu- 
tions remain symmetries for generalized solutions. Moreover, we introduce a gen- 
eralization of the infinitesimal methods of group analysis that allows to compute 
symmetries of linear and nonlinear differential equations containing generalized 
function terms. Thereby, the group generators and group actions may be given by 
generalized functions themselves. 
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1 Introduction 

Symmetry properties of distributions and group invariant distributional solutions 
(in particular: fundamental solutions) to particular types of linear differential 
operators have been studied by Methee (p2|), Tengstrand (p§]), Szmydt and 



Ziemian ([33, 34, pq| , [pq]). A systematic investigation of the transfer of classical 
group analysis of differential equations into a distributional setting is due to Berest 
and Ibragimov (§, ^, [I], ||, fl8|]), again with a view to determining fundamental 
solutions of certain linear partial differential equations. A survey of the lastnamed 
studies including a comprehensive bibliography can be found in the third volume 
of As these approaches use methods from classical distribution theory, their 
range is confined to linear equations and linear transformations of the dependent 
variables. 

Algebras of generalized functions offer the possibility of going beyond these limi- 
tations towards a generalization of group analysis to genuinely nonlinear problems 
involving singular terms, like distributions or discontinuous nonlinearities. In the 
present paper we develop a theory of group analysis of differential equations in 
algebras of generalized functions that allows a satisfactory treatment of such prob- 



lems. This line of research has been initiated in []2§|] and has been taken up in 21]. 
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Applications to different types of algebras of generalized functions can be found in 
@ and 

The plan of the paper is as follows: Section ^ provides a short introduction to 
the theory of algebras of generalized functions in the sense of J.F. Colombeau. In 
section || we consider systems of partial differential equations together with a clas- 
sical symmetry group G that transforms smooth solutions into smooth solutions. 
Assuming polynomial bounds on the action of G, we can extend it to generalized 
functions belonging to Colombeau algebras and ask whether G remains a symme- 



try group for generalized solutions. In section 3A we develop methods based on a 
factorization property of the transformed system of equations. Essentially, poly- 
nomial bounds on the factors suffice to give a positive answer. In the scalar case 
we show this to be automatically satisfied whenever the equation contains at least 
one of the derivatives of the solution as an isolated term. While the conditions of 



section 3.1 concern some mild assumptions on the algebraic structure of the equa- 



tions, section |3.2| develops a topological criterion, applicable to systems of linear 
equations: the existence of a C°°-continuous homogeneous right inverse guarantees 
a positive answer as well. Along the way we give examples of nonlinear symmetry 
transformations of shock and delta wave solutions to linear and nonlinear systems. 
The purpose of section |4| is to develop the general theory, allowing the equations 
and the group action (hence also its generators) to be given by generalized func- 
tions. Using the characterization of Colombeau generalized functions by their gen- 
eralized pointvalues established in |27|] as well as results on Colombeau solutions 
to ODEs, we show that the classical procedure for computing symmetries can be 
literally transferred to the generalized function situation. The defining equations 
are derived as usual, but their solutions are sought in generalized functions. This 
enlarges the reservoir of possible symmetries of classical equations and allows the 
study of symmetries of equations with singular terms. An example is provided by 
a conservation law with discontinuous flux function. 

The remainder of the introduction is devoted to fixing notations and recalling 
some basic definitions from group analysis of differential equations. We basically 
follow the notations and terminology of ]29|| . Thus for the action of a Lie group 
G on some manifold M, assumed to be an open subset of some space X x U of 
independent and dependent variables (with dim(^f ) = p and dim(ZY) = q) we write 
g ■ (x, u) = (E g (x, u), $ g (x, it)). Transformation groups are always supposed to act 
regularly on M. If H 9 does not depend on u, the group action is called projectable. 
Elements of the Lie algebra g of G as well as the corresponding vector fields on 
M will typically be denoted by v and the one-parameter subgroup generated by 
v by i] — > exp(?]v). denotes the n-jet space of M; the n-th prolongation of 

a group action g or vector field v is written as pr( n )g or pr( n )v, respectively. Any 
system S of n-th order differential equations in p dependent and q independent 
variables can be written in the form 

A u (x,u^) = 0, l<u<l. 
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where the map 

A : X x -> R' 
(x,u^) -> (Ai(x, u (n) ),..., Aj(a:, 

will be supposed to be smooth. Hence S is identified with the subvariety 

Sa = {(x,u (n) ) : A(x,u {n) ) =0} 

of A" x «(•»). For any / : Q C AT -» W, 1/ is the graph of / and : = 

{(x, pr( n )/(x)) : x G 0} is the graph of the n-jet of /. 

2 Colombeau algebras 

Already at a very early stage of development of the theory of distributions it 
became clear that it is impossible to embed the space T>'(Q) of distributions over 
some open subset Q of E" into an associative commutative algebra (A(£l),+,o) 
satisfying 

(i) T>'(Q) is linearly embedded into A(Q) and f(x) = 1 is the unity in A(£l). 

(ii) There exist derivation operators di : A(£l) — > A(tt) (i = 1, ... ,n) that are 
linear and satisfy the Leibnitz rule. 

(hi) 9j|x>'(n) i s t ne usual partial derivative (i = 1, . . . , n). 

(iv) °|c(n)xC(n) coincides with the pointwise product of functions. 

(Schwartz's impossibility result, p2[|). Furthermore, replacing C(Q) byC^^(Q) does 
not alter this result. On the other hand, many problems involving differentiation 
and nonlinearities in the presence of singular objects require a method of coping 



with this situation in a consistent manner (cf. e.g. [|2j], |2q1 , ||37| ). By the above, 
the best possible result would consist in constructing an algebra A(Q) satisfying 
(i)-(iii) and 

(iv') °| CO o(Q) xC oc(ft) coincides with the pointwise product of functions. 

The actual construction of algebras enjoying these optimal properties is due to 
J.F. Colombeau (||, H, see also BJ, The basic idea underlying his theory 

(in its simplest - the so-called 'special' - form) is that of embedding the space 
of distributions into a factor algebra of C°°(r2) / (J = (0, 1]) via regularization by 
convolution with a fixed 'mollifier' p G 5(IR' 1 ) with J p(x)dx = 1. In order to 
motivate the definition below let p £ (x) := e~ n p{-) and let u E £'(R n ) (the space 
of compactly supported distributions on R n ). The sequence (u*p £ ) ee j converges to 
u in 2?'(R n ). Taking this sequence as a representative of u we obtain an embedding 
of £>'(R") into the algebra C°°(R) 7 . However, embedding C 00 (R") C P'(R n ) into 
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this algebra via convolution as above will not yield a subalgebra since of course 
(/ * Pe){g * Pe) 7^ (/<?) * Pe i n general. The idea, therefore, is to factor out an ideal 
Af(R n ) such that this difference vanishes in the resulting quotient. In order to 
construct jV(R n ) it is obviously sufficient to find an ideal containing all differences 
(/ * Pe)eei ~ (f)eei- Taylor expansion of / * p £ — f shows that this term will vanish 
faster than any power of e, (uniformly on compact sets, in all derivatives) provided 
we additionally suppose that J p{x)x a dx = for all a£NfJ with \a\ > 1. The set 
of all such sequences is not an ideal in C°°(R ri ) i , so we shall replace C^iW 1 ) 1 by 
the set of moderate sequences Sm^"") whose every derivative is bounded uniformly 
on compact sets by some inverse power of e. 

Thus we define the Colombeau algebra Q{Q) as the quotient algebra £m{&) 
where 

S M (n) := {{u £ ) e£l G C°°(QY : \/K CC ft,Va G N™ 3p G N with 
sup \d a u £ {x)\ = 0{e- p ) as e -» 0} 

xeK 

M{n) := {(m e ) £6 / € r(t2) ; : Vif CC fi,Va G V? G N 
sup \d a u £ (x)\ = 0(e q ) as e -» 0}. 

x<=K 

Equivalence classes of sequences (u £ ) £ ^j in Q(fl) will be denoted by cl[(« £ ) £e /]. 
Q{Q) is a differential algebra containing £'{Q) as a linear subspace via the embed- 
ding l : u — > cl[(u*p £ ) £& i] depending on a mollifier p G 5(M n ) as above, i commutes 
with partial derivatives and coincides with u — > cl[(u) ee /] on V(Q), thus rendering 
it a faithful subalgebra of The functor Q — > £/(f2) is a fine sheaf of differen- 

tial algebras on R™ and there is a unique sheaf morphism i extending the above 
embedding to C°°( . ) <^-> X>'( . ) <^-> C/( . ). t commutes with partial derivatives, and 
its restriction to C°° is a sheaf morphism of algebras. 

We shall also consider the algebra G T (£l) = £ T (Q,)/J\f T (Q) of tempered generalized 
functions, where 

0Af (ft) = {/ G C°°(n) : Va G N™ 3p > sup(l + \x\)~ p \d a f(x)\ < oo} 

£r(fi) = {(u £ ) £eI G (0 M (fi)) 7 : Va G N™ 3p > 

sup(l + \x\)- p \d a u £ (x)\ = 0(e- p ) (e -» 0)} 

= {(u £ ) ££l G {O m {Q)Y : Va G N™ Bp > V 9 > 
sup(l + \x\)- p \d a u £ (x)\ = 0{e q ) (e 0)} 

The map t defined above is a linear embedding of <S'(R n ) into £/ T (R n ) commuting 
with partial derivatives and making 

O c (M. n ) = {/ G C 00 ^™) : 3p > Vq G N™ sup (1 + |x|)- p |d Q /(x)| < oo} 
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a faithful subalgebra. Elements of Om (0) are called slowly increasing. Componen- 
twise insertion of elements of Q into slowly increasing functions yields well defined 
elements of Q. Thus, in Q not only polynomial combinations of distributions (e.g. 
5 2 ) make sense but also expressions like sin(<5) have a well-defined meaning. The 
importance of G T (Q) for our purposes stems from the fact that elements of this 
algebra can even be composed with each other (again by componentwise insertion, 
cf. [jn|], [p0| ), a necessary prerequisite for generalizing symmetry methods, see 
section |4|. Especially in the theory of ODEs in the generalized function context it 
is often useful to consider the algebra Q T (£l x fi') whose elements satisfy (/-bounds 
in the f2-variables and t/ T -bounds in the fJ'-variables (cf. or pQ]). Elements of 
Colombeau algebras are usually denoted by capital letters with the understanding 
that (u £ ) £ £i denotes an arbitrary representative of U G Q. 

Nonlinear operations with distributions in G(£l) depend not only on the distribu- 
tions themselves but also on the regularization procedure used in the embedding 
process. Thus the difference of two representatives (« e ) eg j, {y £ ) £& i of generalized 
functions U resp. V may have D'-limit as e — > without U and V being equal 
in Q(Q). Nevertheless U and V are to be considered 'equal in the sense of distri- 
butions' or associated with each other (U ~ V). Moreover, U is called associated 
with some distribution w if u £ — > w in T>' . If such a w exists (which need not be 
the case, cf. <5 2 ), it is to be seen as the distributional 'shadow' of U. For example, 
all powers of the Heaviside function are associated with each other without being 
equal in the algebra itself. Also, xb = in 2?'(R), so xS ~ in (/(K), but xS ^ in 
Q(M). These examples illustrate a general principle: assigning nonlinear properties 
to elements of the vector space V{Q) amounts to introducing additional informa- 
tion which is reflected in a more rigid concept of equality within £/(Q) compared to 
that in D'(£l). This strict concept of equality allows for much more refined ways of 
infinitesimal modelling. On the D'-level (the level of association) this additional 
information is lost in the limit-process e — > 0. 

Generalized numbers (i.e. the ring of constants in case O is connected) in any of 
the above algebras will be denoted by 1Z. Componentwise insertion of points into 
representatives of generalized functions yields well defined elements of 1Z. 
We note that there exist variants of Colombeau algebras that allow a canonical 
embedding of distributions (indepenent of a fixed mollifier as above). The basic 
idea for constructing these algebras is to replace the index set I by the space of 
all possible modifiers. Our choice of the special variants of Colombeau algebras 
is aimed at notational simplicity. However, all results presented in the sequel 
carry over to the respective full variants of the algebras. Moreover, recently there 
have been introduced global versions of Colombeau algebras, defined intrinsically 
on manifolds and displaying the analogues of (i)-(iv) (with d{ replaced by Lie- 
derivatives with respect to smooth vector fields), sec [Q. For applications of the 
theory to nonlinear PDEs see [24] and the literature cited therein, for applications 
to mathematical physics and numerics, cf. [||, [ 10 1 and |$7fl. 
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3 Transfer of Classical Symmetry Groups 



3.1 Factorization Properties 

The first question to be answered in trying to extend the applicability of classical 
group analysis to generalized solutions concerns permanence properties of classical 
symmetries: Let G be the symmetry group of some system S of PDEs and consider 
S within the framework of Q(Q). Under which conditions do elements of G also 
transform generalized solutions into other generalized solutions? It is the aim of 
this and the following section to answer this question. To begin with, let us fix 
some terminology: 

3.1 Definition Let G be a projectable local group of transformations acting on 
some open set M. C X x U according to g ■ (x,u) = Cz, g (x),& g (x,u)). g is called 
slowly increasing if the map u — > <& g (x, u) is slowly increasing, uniformly for x in 
compact sets, g is strictly slowly increasing if <3> s G Om(M). IfQ^X, U G G{£T) 
and g is (strictly) slowly increasing, the action of g on U is defined as the element 

5 C/:=cl[(($ 9 o( i( ix U£ ))oH- 1 ) £e/ ] (1) 

ofg(E g (n)). 

If U is a smooth function, (|l|) reproduces the classical notion of group action on 
functions. Henceforth we make the tacit assumption that the differential equa- 
tions under consideration are of a form that allows for an insertion of elements of 
Colombeau generalized functions (i.e. the function A representing the equations 
on the prolongation space is slowly increasing). Also, slowly increasing group ac- 
tions are always understood to be projectable. Analogous to the classical setting 
we give the following 

3.2 Definition Let S be some system of differential equations with p variables 
and q unknown functions. A solution of S in Q is an element U G (G(Q)) Q , with 
Q, C X open, which solves the system with equality in (Q(Q)) 1 . A symmetry group 
of S in Q is a local transformation group acting on X xlA such that if U is a 
solution of the system in Q, g G G and g-U is defined, then also g-U is a solution 
of S in Q . 

Let us take a look at the transition problem from classical to generalized symmetry 
groups on the level of representatives. Thus, let G be a slowly increasing symmetry 
group of some differential equation 

A(x,n( n )) = 0. (2) 

This means that if / is a classical solution, i.e. if A(x, pr^ n ^/(x)) = for all x then 
also A(:r, pr( n )(g ■ f){x)) = 0. Now let U G G{Q) be a generalized solution to (||). 
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Then for any representative (u e ) eg / of U there exists some (re e ) ee j € jV(fi) such 
that for all x and all e we have 

A(x,pr (n) u £ (x)) =n e (x). (3) 

In particular, the differential equation (||) need not be satisfied for even one single 
value of e. This basic observation displays quite fundamental obstacles to a di- 
rect utilization of the classical symmetry group properties of G in order to obtain 
statements on the status of G in the Colombeau-setting. Therefore we have to 
derive properties of symmetry groups that are better suited to allow a transfer to 
differential algebras. The starting point for our considerations is a slight modifica- 
tion of a well known factorization property of smooth maps (cf. |2J| , Proposition 
2.10): 

3.3 Proposition Let F be a smooth mapping from some manifold M to M. k (k < 
n = dim(M)), let f : (—n ,rj ) x M — > R be smooth and suppose that f(rj, .) 
vanishes on the zero set Sf of F, identically in n. If F is of maximal rank (= k) 
on Sf then there exist smooth functions Q\, . . . , Q k '■ {—f]oi Vo) x M — * R such that 

f(n,m) = Q 1 (r),m)F 1 (m) + . . . + Q k (r],rn)F k (m) 

for all (j],m) € (— f?o,f?o) x M. □ 



We are mainly interested in the following application of Proposition 3.3: 
3.4 Theorem Let 

A u (x,u {n) ) = 0, \<v<l (4) 

be a nondegenerate system of PDEs. Let G = {g^ : n € (—r] ,rj )} be a one 
parameter symmetry group of (0j and set g v ■ (x,u) = (E v (x,u),^jj(x,u)). Then 
there exist C°° -functions : (—r] ,r] ) x V — > R (1 < fx, v < I, V an open subset 
of M {n) ) such that if u : n C MP -> R<? is smooth and g v u exists we have 

A v (E v (x, u(x)),p^ n \g v u)(E v (x, u(x)))) = 
i 

= Yl Q At i,(??,a;,pr (n) ii(x))A^(x,pr( n )it(x)) (5) 
on the domain of g v u for 1 < v < I. 

Proof. Denote by z the coordinates on M^ n \ That g„ is an element of the 
symmetry group of the system is equivalent with 

A(z) = o => K{w (n) g n {z)) = o {l<u<l) 

for all rj and z such that this is defined. A is of maximal rank because (Q) is nonde- 
generate. Hence, by Proposition 3J3 there exist C°°-functions : (—r) , rj ) xV ^ 
R (1 <//</, V an open subset of M^) such that 

I 

A„(prM<7„(*)) = E Qrfa z )^ z )- ( 6 ) 
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Now for a smooth function u : f2 C R p — > R q as in our assumption and x G $1 we 
set 

z u (s) := (x.prWufc)) G ^W. (7) 

Then by definition pr( ra ) g r] (z u (x)) = (E v (x,u(x)),pr ( - n \g v u)(E rj (x,u(x)))), so the 
result follows. □ 

For a single PDE A(x,pr^u) = 0, equation (||) takes the simpler form 

A(E v (x, u(x)), pr (n) (^u)(H, ? (x, = Q(r/, x, pr (n) n(x))A(a;, pr (n) u(x)). 

(8) 

Theorem |3.4| will be one of our main tools in transferring classical symmetry groups 
of (systems of) PDEs into the setting of algebras of generalized functions. 

3.5 Proposition Let r\ — > g„ be a slowly increasing one parameter symmetry 
group of IfP^ := (Q^(r?, H_ r ,( . ),pi^u £ (E^, ri ( . )))) ee j belongs to £ M {Q) for 
1 < A*, v < ^ a^rf every (u e ) ee / G £m(^)> ^ien rj ^ g v is a symmetry group of ^) 
in Q as well. This condition is satisfied if 

is slowly increasing in the -variables, uniformly in x on compact sets for 1 < 
fi, v < I and every r). 

Proof. It suffices to observe that (|5|) gives 

A v {x,pr^{g v u){x)) = 
i 

= Y, Q»»(V, 3-„(s), prW u (H_ r ,(x)))A M (3_ 7? (x), pr^H^ (*))). 

For any solution U G G(Q) with representative u = (u e ) £< zj, this expression is in 
A/"(0) since G for each fi, v, and every A /J (H_ J? ( . ), pr^ n )u(H_ J? ( . ))) is 

in A/"(0) because U is a solution and is a diffeomorphism. □ 

3.6 Example The system 

U t + UU x = 

V t + UV x = (9) 

u\ {t=0} =u o , v\ {t=0} =v o 

may serve as a simplified model for a one-dimensional, elastic material of high 
density in a nearly plastic state. It was analyzed in [25], where solutions U,V G 
£?s, s (R x [0, oo)), J7 , V G ^ S)9 (M) were constructed and studied (G s ,g is a variant of 
the Colombeau algebra with global instead of local bounds). In the following we 
present some applications of the above results to this system (for a more detailed 
study, see 0). For U' a > ® has a unique solution (^V) in </ s , g (M x [0,oo)) 
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with d x U > 0. We consider solutions in Q Stg (M. x [0, oo)) with initial data U Q (x) = 
u L + (u R — u L )H(x) and V (x) = v L + (v R — v L )H(x), where H is a generalized 
Heaviside function with H' > 0, i.e. H is a member of G St9 (M.) with a representative 
(h e ) e< zj coinciding with the classical Heaviside function Y off the interval [— e, e]. 
For u L < u R the solution (U, V) is associated with the rarefaction wave 



u(x, t) 



U L , X < U L t 

I , u L t <x < u R t (10) 
Ur, U R t < X 



v(x, t) 



v L , x < U L t 

2B=^) f + u L t<x<u R t (11) 

V R , U R t < X 



However, choosing different generalized Heaviside functions for modelling the ini- 
tial data U , respectively V Q we may obtain a superposition of the rarefaction wave 
( |l0|) in u with a shock wave 

v{x,t) = v L + (v R -v L )Y{x-ct) (12) 

with arbitrary shock speed c, u L < c < u R . We are going to construct a one 
parameter symmetry rj — > g„ of (||) which transforms any of the solutions ( |Tll) , ([12] ) 
into a shock wave solution as rj — > ±oo. For this we employ the two-dimensional 
Lorentz-transformation (77, (x,t)) — > (xcosh(r/) — tsinh(rj), — xsmh(rj) + tcosh(r/)) 
with infinitesimal generator X Q = —td x —xdt- Then X := X + (u 2 — l)d u generates 
a projectable one-parameter symmetry group of (||). Assuming that — 1 < u L < 
n R < 1, we can extend the solution (U,V) to the region Q = M 2 \ {(x,t) : t < 
0, u R t < x < u L t} by the method of characteristics applied to representatives. 
Then the Lorentz-transformed solutions 

u £ (x,t) = — tanh(?7 — Artanh(u £ (x cosh(r/) + t sinh(?7), 

x sinh(r?) + t cosh(?7)))) (13) 
v £ (x,t) = v £ (x cosh(rj) + tsinh(ry),rcsinh(r/) + tcosh(r/)) (14) 

(with Artanh the inverse of tanh) are well defined at least on K x (0, 00). The 
factorization property (||) in this case reads 

(d t u £ + u £ d x u £ )(x,t) = (15) 
{(dtu £ + u £ d x u £ ) / (cosh 3 (Art&nh(u £ — 77)) cosh(Artanh(-u e )))) (E~ 1 (x, i)) 

and similarly for the second line in @, demonstrating that (U ,V) is again a 
solution. For each rj, U is associated with a piecewise smooth function which 
converges to =Fl as rj — > ±00. Observing that the coordinate transformations in 
(|l3|), (|i~4|) approach boosts in the directions (=Fl, 1) as rj — > ±00, we see that the 
functions associated with V converge to the shock wave t> L + (v R — v L )Y(x db i) as 
77 — > ±00, for whatever solution V given in ( JTT| ) or (p^). 
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Although Proposition 3.5 provides a manageable algorithm to determine if classical 
symmetry groups carry over to generalized solutions it would certainly be prefer- 
able to have criteria at hand that allow to judge directly from the given PDE if the 
factors P^ u behave nicely (given slowly increasing group actions). The first step 
in this direction is gaining control over the behaviour of the map z — » j>A n ^ g v (z), 
defined on M^ n \ 

3.7 Proposition If rj — > g v is a (strictly) slowly increasing group action on Ai 
then z — > pi^gri(z) is (strictly) slowly increasing as well. 

Proof. Let N := dim(.M (n) ). For z = (zi, . . . , z p , z p+1 , . . . , z g , . . . , z N ) e M {n) 
we choose some smooth function h : X — » U satisfying z = Zh(z%, . . . , z p ), with 
Zh(x) as in (Q). Then we set x := (z±, . . . , z p ), u = (z p+ i, . . . , z q ), x = r, v (x) 
and u = <&r)(x, u). By the definition of prolongued group actions we have to find 
estimates for every 

A a :=((*„ o(idxh))o ~_„) (s) (i) (16) 

(where (s) denotes the derivative of order s) in terms of z. The above formula 
contains the components of pA n ^g(z) of order s (s < re). Note that the particular 
choice of h has no influence on (jig), i.e. A s depends exclusively on z. To compute 
A s explicitly we use the formula for higher order derivatives of composite functions 
(see [11 1). Denoting by T m the group of permutations of {1, ... , m} we have: 

A s (r u ...,r s ) = J2J2J2liu^°(idx h)f\(x))(t u . . . , t % ), (17) 



where 



h = H^ ) (a;)(r ff ( 1) ,... ,r CT(A;i )), . . . ,U = E™(a?)(r CT ( s _ fei+ i), • • • ,r a (s))- 



and 



(((#, o (id x (x)(t 1 ,...,i l ) = EE E ^fo • • • > o> 

|iM 

(18) 

where 

si = (idx h)^(x)(t T{1) ,...,t T{ll) ), ..., Sj = (idx h)^(x)(t Tii _ lj+1] ,...,t T(i} ). 

Each s m consists of sums of products of certain t T ^) with certain zi and an anal- 
ogous assertion holds for the $fj (x, u)(s\, ... ,Sj). Hence from (|l7|) and fll8|) the 
result follows. □ 

Returning to our original task of finding a priori estimates for the factors P^ v , even 
with the aid of Proposition 3.7 we still need some information about the explicit 
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form of the to go on. In general this seems quite difficult to achieve. However, 
there is a large and important class of PDEs that allow a priori statements on the 
concrete form of the factorization. Namely, we are going to show that each scalar 
PDE in which at least u or one of its derivatives appears as a single term with 
constant coefficient belongs to this class. 

Consider a scalar PDE A(x,u^) = together with a symmetry group r\ —* g v . 
Then we have 

A(z) = => A(pr^g v (z)) = 

Set F(z) := A(z), f(z) := A(pr^g v {z)) and N = dim(M^). Suppose that in 
a neighborhood of some z with F(z) = we have j£- > for some 1 < k < N. 
Then by the implicit function theorem, locally there exists a smooth function 
tp : R N ~ 1 — > K such that in a suitable neighborhood of z we have 

F(z) =0 e> z k = iP(z'), 

where z' = {z\, . . . , z k , . . . , zn) (meaning that the component z k is missing from 
z'). It follows that 

F(z) = (z k - ip{z)) / — (zi, . . . ,z k -i,Tz k + (1 - r)ip(z'), . . . ,z N )dr, 



o dz k 



and on the other hand 



r 1 df , 
f(z) = {z k -ip{z')) J — (z 1 , . . . ,z k -i,Tz k + (1 - r)ip(z'), . . . ,z N ) dr. 

Thus in the said neighborhood we have 

tt ^ t?i J° ^(^■■■^k-i,rz k + (l-r)iP(z'),...,z N )dr 

j(z) = F(z) 1 " (19) 

Jo §k( Zl >' ' ' ' Zk ~ 1 ' TZk + ^ ~ T )^ Z ')^ • • • ' z N)dT 

provided the denominator of this expression is 7^ 0. In particular, if for some 
constant c/Owe have Jj- = c in a neighborhood of z then (|l9| ) simplifies to 



f(z) = -F{z) 1 ° J 
c 



/ — (zi,..., z k -i,Tz k + (1-t)i/)(z ),..., ZN)dr (20) 
After these preparations we can state 



3.8 Theorem Let rj — > g v be a slowly increasing symmetry group of the equation 
A(x,u^ n ') = 0. Set N = dim{M^ n >) and suppose that j^- = c 7^ for some 
p + 1 < k < N. Then n — > g v is a symmetry group of A(x, u^) =0 in Q. 

Proof. Without loss of generality we may assume c = 1. Using the above nota- 
tions we have F(z) = z k — ip(z'), so (P0|) implies 



f df 

f(z) = F(z) J -Q^-izii ■ ■ ■ , z k -i,Tz k + (l-T)(z k -F(z)), . . . ,z N )dr =: F(z)Q(rj,z) 
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From Proposition 3.7 we know that z — * f(z) is slowly increasing in the u^ n '- 



variables (i.e. in those z% with i > p), uniformly in x = (zi,... ,z p ) on compact 
sets. Since F is slowly increasing we infer that Q(r], z u (x)) G £m{&) f° r an Y 
u G £m{&) (with z u as in (0)). Finally, 

A(x, p^ n Hg vU )(x)) = A(E-r,(x),pr^u(E^(x)))Q(ri, E^(x),pr^u(Z_ v (x))). 

Since E_ v is a diffeomorphism, it follows that if U = cl[u] solves the equation, so 
does grjU. □ 

As the proof shows, we can drop the assumption p+1 < k if we require the group 
action to be strictly slowly increasing. It is clear that many PDEs satisfy the 
requirements of Theorem |3.8[ For example, in the Hopf equation A(x, t, u, u Xl u t ) = 
ut + uu x or A(zi, . . . , 25) = 2:5 + Z3Z4 one can take k = 5. Note however that 
not every symmetry group of this equation is automatically slowly increasing. 



Theorem 3.5 constitutes a useful tool for transferring classical symmetry groups 



to Colombeau algebras. 

3.9 Example We consider the initial value problem for the nonlinear transport 
equation 

U t + A • V X U = f(U) 

u\ {t=0} =u o 

with t G R,x,X G R™. It has unique solutions in £(R n+1 ), given U Q G G(R n ), 
provided / G Om is globally Lipschitz (see [p^l ). If in addition / is bounded and the 
initial data are distributions with discrete support, say Uq(x) = ^ • a^^W [x — ^) 
with £j G R n ,i G Ng, then the generalized solution is associated with the delta 
wave v + w where 

v {x,t) = ^a ij 5 {i) {x-\t-i j ) (22) 

hj 

and w is the smooth solution to wt + A • V x w = f(w), w(0) = 0. 
The vector field X = cf(u)d u generates an infinitesimal symmetry of (^l|) for arbi- 
trary c G M. With F(u) := J du/f(u), the corresponding Lie point transformation 
is 

(x,t,u) -» (x,t,u) = (x^^F- 1 (crj + F(u))). (23) 

This provides a well-defined nonlinear transformation of the generalized solution 
U G Q(R n+1 ), provided that the right hand side in (Ell) is slowly increasing. 
In the example 

U t + X- V X U = tanh(i7) (24) 

the generalized solution is associated with v(x,t) and w vanishes identically. Ap- 
plying (]23| ) we obtain (due to Theorem |3.8| ) the new generalized solution 

U[x, t) = Arsinh (e cr? sinh(?7(x, t))) (25) 
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(with Arsinh the inverse of sinh). We are going to show that U is still as- 
sociated with the delta wave v in (|22]), To simplify the argument we assume 
n = 1, A = and Uq(x) = 5^\x). Representatives of U resp. U are u £ (x,t) = 
Arsinh(e*sinh( y O£*' ) (x))) and u £ (x,t) = Arsinh(e c,?+ *sinh(p^(x))). For ip G V(R 2 ) 
we have 

P £ := f f u £ (x,t)ifj(x,t)dxdt = 
= fff ( ! l 8(e c i+ t ,ae- l - 1 p( i \x))dae~ t p ( - l \x)ip(ex,t)dxdt 

where 9(a,y) := ^Arsinh(asinh(y)) for a > 0, y € R. Since 9 is bounded by 

max(l,a) and lim^i^^ 9(a,y) = 1 it follows that I® — > J ip(0,t)dt, so U is associ- 
ated with the delta function on the i-axis, as desired. For i > 1 we write 

4 = I J 'f<5( (, eCV+1 \oe~ l - A p^{x)) - l)dae- i p^{x)ip{ex,t)dxdt+ 
+(-iy J J p{x)d l x ip(ex,t)dxdt 

Here the second term converges to (—1)* J 9^(0, t) and the first term goes to zero 
since Jq \9(a, ay) — l\da < i 1 ^ (1 — e~^) for y ^ 0. This proves the claim for 
p € T>(M.). For p € <S(R) splitting the x-integral into one from — to 4g and one 
over I a; I > -4= gives the same result. 

3.2 Continuity Properties 

In this section we work out a different strategy for transferring classical point sym- 
metries into the ^-setting. This approach, suggested in 28], consists in a more 



topological way of looking at the transfer problem by using continuity properties 
of differential operators. As we have pointed out in the discussion following (||), 
the main obstacle against directly applying classical symmetry groups componen- 
twise to representatives of generalized solutions is that the differential equations 
need not be satisfied componentwise. However, there are certain classes of partial 
differential operators that do allow such a direct application. Consider a linear 
partial differential operator P giving rise to an equation 

PU = (26) 



in Q and let G be a classical slowly increasing symmetry group of (26). Further- 
more, suppose that P possesses a continuous homogeneous (but not necessarily 
linear) right inverse Q. If U = cl[u] is a solution to (|26l ) in G(fl) then there exists 
some n € Af(£l) such that 

Pu = n. 

Since Q is a right inverse of P this implies 

P{u £ - Qn £ ) = Ve € /. (27) 
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Also, Qn 6 AA(f2) due to the continuity and homogeneity assumption on Q. If 
g G G, (|^) implies 

P( 5 (n e - Qn s )) = VeGJ. 

By definition, 

P(gU) = cl[P(gu)} = cl[P(g(u - Qn))}, 

so gU is a solution as well. Summing up, G is a symmetry group in Q. The 
following result will serve to secure the existence of a right inverse as above for a 
large class of linear differential operators. 

3.10 Proposition Let E, F be Frechet spaces and A a continuous linear map 
from E onto F . Then A has a continuous homogeneous right inverse B : F — > E. 



Proof. See j23|, p. 364. □ 
From these preparations we conclude 

3.11 Theorem Let 

A u (x,u (n) ) = 0, u = l,...,l 

be a system of linear PDEs with slowly increasing A u and let n — > g v be a slowly 
increasing symmetry group of this system. Assume that the operator defined by 
the left hand side is surjective (C°°(0))' — > (C°°(£l)) 1 . Then n — > g v is a symmetry 
group for the system in Q{£1) as well. □ 

The assumptions of Theorem |3.11 are automatically satisfied for any linear partial 
differential operator with constant coefficients on an arbitrary convex open domain 
(see 0, 10.6). 

3.12 Example The system of one-dimensional linear acoustics 

Pt + U x = , , 

u t + p x = o. (28) 



is transformed via U = V — W,P = V + W into 

(29) 



v t + V x = 



w t -w x = 0. 

Using the infinitesimal generators &(v)d v + *$>(w)d w ($, ^ arbitrary smooth func- 
tions) of (|29| ) we obtain symmetry transformations for ( ^8|) of the form 



U = F- 1 (r? + F{\(P + [/))) - G- 1 (e + G(\(P - U)f) 
P = F- 1 L + F{ l -{P + U))\ + G- 1 (e + G( l -{P - U)) 

with arbitrary diffeomorphisms F, G. Since (p§|) satisfies the assumptions of The- 
orem 3.11 on £2 = R 2 it follows that any slowly increasing transformation of this 
form is a symmetry of (p8|). In particular, this includes nonlinear transformations 
of distributional solutions, cf. Example 3.13; . 
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In the remainder of this section we discuss the interplay between symmetry groups 
and solutions of PDEs in the sense of association. Consider 



A u (x,u {n) ) « 0, l<u<l (30) 

in Q. A slowly increasing symmetry group of the corresponding system 

A{x,u {n) ) = 0, \<v<l 

is called a symmetry group in the sense of association if it transforms solutions of 
( |30| ) into other such solutions. The first question to be answered in this context 
is whether one can derive conditions on the form of the factorization (||) that will 
yield symmetry groups in the sense of association. It is clear that a sufficient 
condition is to suppose that Q depends exclusively on rj and x. Distributional 
solutions to linear PDEs arise as a special case of ( |30[ ) and have been treated 
in There, the validity of equation (|8|) with Q depending on r\ and x only 
is actually used to define symmetry groups in T>' . In order to remain within 
the classical distributional framework, the admissible group transformations in ||] 
are restricted to projectable ones acting linearly in the dependent variables. On 
the other hand, the method developed there is even applicable to linear equations 
containing distributional terms which allows to use invariance methods to compute 
fundamental solutions. 

Second, if u is a solution to A(x,u^) = in possessing an associated distri- 
bution, one may ask for which group actions g this implies that gu as well possesses 
an associated distribution. This is certainly the case for admissible transforma- 
tions in the above sense. On the other hand, we have already seen in Example |3.9| 
that even genuinely nonlinear symmetry transformations may preserve association 
properties. 

The next example shows that nonlinear group actions may transform distributional 
solutions in Examples 3Jj and |3.12| into more complicated distributional solutions 



or into generalized solutions in Q(M. ) not admitting associated distributions. 



3.13 Example We consider the equation Ut + XU X = arising in (|2l| ) with n = 1 
or in (29). We have already observed that U = F~ 1 (rj + F(U)) defines a symmetry 



transformation for arbitrary diffeomorphisms F. Here we take F G C°°(R), F' > 0, 
F(y) = sign(y)^/jy| for \y\ > 1. We wish to compute U when U G Q(M. 2 ) is a delta 

wave solution U(x,t) w 5^(x — \t). We take U as the class of p^\x — At) with 
p G T>([—1, 1]). We have when rj > 0: 

(i) If i = 0, that is U rj S(x - Xt), then U rj F~ l {r] + F(0)) + S(x - At); 

(ii) If % = 1, that is U rj S'(x - At), then 

U w F-^ri + F(0)) +2 V J y/\/7(y)\dy 5(x - At) + 8'{x - At); 

(iii) If i > 2 then U does not admit an associated distribution. 
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To see this, we may assume that A = and write a e (x) := p £ ■ {%) for brevity. Note 
that F~ l (y) = sign(y)y 2 for \y\ > 1. Let A £ = {x £ [—e,e] : |a e (x)| > (rj + l) 2 }. 
If x 6 A £ and a £ (x) > then rj + F(a £ (x)) > 1 and F -1 ^ + F(a £ (x))) = rj 2 + 
2r]^Ja £ ~{x) + a e (x). Also, if x 6 A e and a e (x) < then rj + F(a e (x)) < —1 and 
F~ x {r] + F{a £ (x))) = -rj 2 + 2r/ v / |a £ (x)| + a e (x). The functions F _1 (ry + F(a e )), 
|a e (x)| and ^/|a e (a;)| are bounded on the complement of A £ . Thus 

/ j £ _ £ F- 1 ^ + F(a £ (x)))Tp{x,t)dxdt = 

= J J A£ (±7] 2 + 2riy/\a e (x)\ + a s (x))1)(x,t)dxdt + 0(e) = 

= J j £ _ e {2r,^/\a £ {x)\ +ae(x))il>(x,t)dxdt + 0(e) 

while 

/ / F- 1 (rj + F(a £ (x)))ip(x,t)dxdt -» F^ + i^O)) / [ ip(x,t) dxdt 

J J\x\> £ J J 

It follows that i ?_1 (r/ + i ? (a £ (a;))) converges in P'(1R 2 ) if and only if 2r/y / |a e | + a e 
admits an associated distribution. A simple computation yields the particular 
results (i), (ii), (hi). 



4 Generalized Group Actions 

Although the methods introduced in the previous sections enable an application 
of large classes of classical symmetry groups to elements of Colombeau algebras, 
they are but the first step in a theory of generalized group analysis of differential 
equations. In this section we develop an extension of the methods of group analysis 
that will allow to consider symmetry groups of differential equations whose actions 
are generalized functions themselves. 



4.1 Generalized Transformation Groups 

Simple examples indicate the necessity of extending the methods of group analysis 
of PDEs to equations involving generalized functions themselves: 

4.1 Example Considering ([2~I|) in Q T with a generalized function / = cl[(/ e ) ee /] 
£ Q T we can apply the classical algorithm for calculating symmetry groups com- 
ponentwise to the equations 

d t u £ + A • V x u E = f £ (u £ ) 

thereby obtaining infinitesimal generators with generalized coefficient functions. 
Thus the question arises in which sense such generators induce symmetries of the 
differential equation. More generally, one can consider differential equations in Q T 
of the form 

P( X) [/M) = 
where P is a generalized function. 
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As is indicated by Example 4.1, composition of generalized functions will inevitably 
occur in a generalization of group analysis. For this purpose, we shall apply suitable 
variants of Colombeau algebras for the following considerations, namely Q T (M. n ) 
and G T (R x E n ) = § T (R 1+n ). 

4.2 Definition A generalized group action on M. n is an element $ of (£/ r (^ 1+n )) n 
such that: 

(i) $(0, .) =idin {G T (R n )) n . 

(ii) $(771+773, •) = Hvi,HV2, •)) m (£ T (K 2 +")r. 

Before we turn to an infinitesimal description of generalized group actions let us 
shortly recall some basic definitions from [27| that are needed for a pointvalue 



characterization of generalized functions which in turn plays a fundamental role 
in the following considerations. Thus for any open set Q C R n we set 

n M ■■= {(x £ )eei eQ 1 :3p> 3 7]>0 \x £ \ < (0 < e < ??)}. 

On £Im we define an equivalence relation by 

[x £ )ea ~ (jfcW & v 9 > 3 V > \x £ - y £ \ < e q (0 < e < 7?) 

and set Q := Om/ ~- & is called the set of generalized points corresponding to Q,. 
The set of compactly supported points is defined as 

Q c = {x € ft : 3 representative (x £ ) £g / 3K CC £1 3rj > : x £ £ K, e G (0,7/)}. 

Note that for O = K we have O = 7£. Theorems 2.4, 2.7 and 2.10 of |27|] establish 



that elements of £7(f2), £7 T (f2) or £/ T (f2 x O') are uniquely determined by their 
pointvalues in O c , f2, or f2 c x fi', respectively. For the theory of ODEs in the 
Colombeau framework we refer to 

4.3 Definition Let £ = ■ ■ ■ ,£ n ) € (£/ r (K n )) n . T/ie generalized vector field 

n 

X = ^2 ii{x)d Xi is called Q -complete if the initial value problem 

8=1 

X (t) = e(x(t)) 

is uniquely solvable in Q(R) n for any x a £ lZ n and any t Q € M. 



4.4 Definition Let & be a generalized group action on R n and set 

$(77, .) e(G T (R n )) n . 



^ dr\ 



o 

n 



7/ the generalized vector field X = £,i(x)d Xi is Q-complete, then X is called the 

i=l 

infinitesimal generator of In this case, $ is also called Q-complete. 
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By [16], every generalized vector field with (^-components whose gradient is of 
L°°-log-type is (/-complete. The notion of infinitesimal generator is well-defined 
due to 

4.5 Proposition Every Q-complete generalized group action is uniquely deter- 
mined by its infinitesimal generator. 

Proof. Let 3>', $" be two (/-complete generalized group actions with the same 
infinitesimal generator X = Y^i=i ^i( x )d Xi - Then both functions satisfy 

-^-$(77,0?) = — $(rj + (J,,x) = — $(/jL,$(r), x)) = £($(77, x)). 
dr) dfx d/i o 

Now given any x G 1Z U , it follows that both 77 — > $'(77, x) and 77 — > $"(77,5?) solve 
the initial value problem 

x(rf) = Z(x(rf)) 
x(0) = x 

By assumption this entails that $'( . ,x) = <J>"( . , x) in ((/(R)) n . Consequently, 

$' (77, x) = $ ' (77, x) 

for all 77 G 7£ c and all x G 7£ n . The claim now follows from |2^|, Theorem 2.10. □ 



As in the classical theory, we are first going to investigate symmetry groups of 
algebraic equations: 

4.6 Definition Let F G (/ T (K n ) and let $ be a generalized group action on M. n . 
$ is called a symmetry group of the equation 

F(x) = 

in g T (R n ) if for any x G lZ n with F(x) = G K it follows that 77 -» F($(t/,x)) = 
in Q(M) (or, equivalently, F(&(rj,x)) = in TZ for every 77 G 1Z C ). 

A characterization of symmetry groups of (generalized) algebraic equations in 
terms of infinitesimal generators is provided by 

4.7 Theorem Let F G Q T (R n ) be of the form 

F\X\ , • • • , X n ) — Xj f\X\y . . . , Xj_i , Xj-|_l , . . . , x n ) 

/or some 1 < i < n and f G ^ T (M n_1 ). Lei & be a Q-complete generalized group 
action with infinitesimal generator X = Y2i=i &( x )® x i an ^ su PP ose that x' — > 
£(x',/(x')) defines a generalized vector field on R n_1 stic/i £/iai i/ie corresponding 
system of ODEs possesses a flow in ((/ r (R 1+ ( n_1 ))) n_1 . The following conditions 
are equivalent: 



18 



(i) $ is a symmetry group of F(x) = 0. 

(ii) Ifx£TZ n with F(x) = G ft it follows that X(F)(x) = in ft. 

Proof, (i) =>• (ii): Consider the function (77, x) — > F(<f>(r),x)) G Q T (M. 1+n ). We 
have 

d n dF 

— F($(77, x)) = ^ ^r(*(»7, *)) = ^X^, *)), 



so that J-| o F($(r/,x)) = X(F)(x) in £ T (M n ). Let x S F such that F(x) = 0. 
Then F($(.,x)) = in Q(M). Thus ^| o F($(r/, z)) = in ft which means that 
X(F)(x)=0mK. 

(ii) (i): We assume F(x\, . . . ,x n ) = x n — f{x\, . . . ,x n -i) and abbreviate (xi, 
. . ., x n -i) by x' . Our first claim is that 

n-1 

= ^^(x'Jix'^djfix') in GAR^ 1 ) 
i=i 

Indeed, if x' € ft™" 1 then F(x',f(x')) = in ft. Hence X(F)(x', /(x')) = in 
7?. for all x' by our assumption. Our claim now follows from |27]], Theorem 2.7. 
Consider the following system of ODEs in Q T : 

x 3 (t) = tjix'Jix')) (j = l,...,n-l) 
x'(0) = a' E ft n_1 

By our assumption, this system has a flow (rj,a') — > (hi(rj,a'), h n -i(r],a')) 
in (^ r (R 1+ ( n-1 ))) n_1 . Set ^(77, a) := /(/ii(t?, a'),. . . ,/in-ifa, a'))- Then 5n (0,a) = 
/(a') and 

fl^a) = (51(77,0), ■■ ■ ,9n(r],a)) := (hi(rj, a'), . . . , K-l(v, a'),9n(v, a )) 

is in (G T (M. 1+n )) n . If a £ ft™ then F(g(rj,a)) = in ft for all 77 € ft c . Therefore, if 
we can show that g( . ,a) = $( . , a) in (£/(R)) n for all a with -F(a) = 0, the proof is 
completed. Now we have gj (77, a) = £j {g\ (77, a) , . . . ,g n (r],a)) for 1 < j < n — 1 and 

g n (rj,a)= J2 9^(01(77, a), . . . , g n -i(v, a))9i(v, a ) = 

i=l 

= Cn(5i(7?, a),..., / (51(77, a), . . . , 5^-1(77, a))) = ^(0(77, «))• 

If F(a) = in ft then a n = /(a'), so that 5(0, a) = (a', /(a')) = = *(0, a). Thus 
g(. ,a) and <&( . , a) solve the same initial value problem. Since X is (/-complete, 
the claim follows. □ 
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4.2 Symmetries of Differential Equations 



In this section we are going to apply the above results to symmetry groups of 
differential equations involving generalized functions. To this end, we will first 
have to define generalized group actions on generalized functions. Once we have 
done this, by a symmetry group of a differential equation we will again mean a 
group action that transforms solutions into other solutions. Thus, from now on 
we will exclusively consider group actions on some space MP x M. q of independent 
and dependent variables. 



4.8 Definition A generalized group action $ G (G T (J 
jectable if it is of the form 



x W+i))p+Q is called pro- 



$(77, (x, u)) = (E v (x), V v (x, u)), 
where E G (<? T (K x RP)J> and V G (C? r (R x RP+«))«. 
The group properties in this case read: 



in g T (RP) Vt?!,^ g n c . 

ty m+V2 (x,u) = ^ Vl (E m {x)^ m {x,u)) in g T (M. p+q ) Vtji,tj2 G K c . 



In particular, we have 



id in G T (R P ) V?? G TZ C 



(31) 
(32) 



(33) 



An adaptation of Lie group analysis to spaces of distributions faces the fundamen- 
tal problem that while the methods of classical Lie group analysis of differential 
equations are geometric in the sense that group action on functions is defined via 
graphs, in classical distribution theory there is no means of defining graphs of dis- 



tributions. However, due to the pointvalue characterization obtained in [27| this 
problem can be dealt with in a satisfactory manner within Colombeau algebras: 

4.9 Definition Let U G {Q{RP)) q and V G {g r {W)) q . The graphs of U and V 
are defined as 

Tu := {(x, 17(5)) : x G K p c } 
T v := {(x,U(x)) : x G 1Z P }. 



It follows directly from [27|, Theorems 2.4 and 2.7 that any generalized function 
is uniquely determined by its graph. Our next aim is to define generalized group 
actions on generalized functions. As in the classical case this is done geometri- 
cally, i.e. by transformation of graphs. The following result is immediate from the 
definitions: 
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4.10 Proposition Let U G (G T (M?)) q and letQ be a projectable generalized group 
action on W x R q . Then $> v (Tu) = T<s> v (u) ^ n Ti p+q for each r\, where ® V (U) 
denotes the element 

x -> ^ v (E_ v (x), U o E- V (x)) 

f(g T (w)y. □ 

We are now able to give a geometric characterization of solutions of PDEs in Q T . 

4.11 Proposition Consider the system of PDEs 

A u (x,U {n) ) = l<v<l (34) 
in G T {R?)) q (where A G (£ r ((R p x R q )^)) 1 ). Set 

5 A := {z G 7e (n) : A„(i) = (1 < ^ < /)}. 
T/ien C7 G (£? r (R p )) 9 is a solution of the system iffT-p T ( t „) U C 5a- 



Proof. This follows immediately from [27], Theorem 2.7. □ 

Prolongation of generalized group actions can be handled in a similar fashion as 
in the classical theory. Thus, let $ be a projectable generalized group action 
on R p x R q . We want to define the re-th prolongation pr( n )<& as a projectable 
generalized group action on (R p x R9)W. Let z G (R p x R9) (n) and choose ft G 
O m (K p ) 9 such that Ox,.. . , z p , pr( n )/i(zi, . . . ,z p )) = 2. Now set 

prto$fa,z) := (^(z!,...^^,?^)^^))^^,...,^))). (35) 

Using for /i a suitable Taylor polynomial, it follows that pr^$ G (<5 T (R x (R p x 
R 9)(n))JV ( w here TV = dim((RP+9)( n ))). Mor eover, the definition does not depend 
on the particular choice of h, which follows exactly as in the classical case. 

4.12 Lemma Let z G (Jl p x TZ q )^ and assume that U G (G T (R p )) q satisfies (ii, 
. . . , z p , px^ n 'U(zi, . . . , Zp)) = z. Then 

prW*(»7,2) = . . . ,2p),pr {n \%(U) )(E,(z u . . . ,z p ))) Vr? G ft c . (36) 

Proof. Let U = cl[(u e ) eg /] and choose a representative (^ e ) ee j of if such that 

(Zu, . . . , Z p£ ,pV^U £ (z l£ , . . . , Z pe )) = Z £ Ve. 

Using the chain rule as in Proposition |3.7|, it follows that the right hand sides of 



fl35j) (with z replaced by z) and of (36) have the same representative (depending 
exclusively on (z e ) e6 /). □ 

4.13 Proposition pr( n )<J> is a generalized group action on (RP x «<?)("). 
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Proof. Property |4.2| (i) is clearly satisfied. Concerning (ii), according to 27], 
Theorem 2.7 it suffices to show that 

pr W$( 7?1 + m ,T) =pA n ^( m ,pi ( - n ^(r ]2 ,I)) Vr/1, 772 € K c , We (K p x K«)W. 

Choose some U G (^ r (R p ))' ? with . . . , z p , pA n W(zi, . . . ,z p )) = z. Then due to 
Lemma 4.12 we have 



By 



pr(")$(r ?2 ,i) = (H 7 , 2 (? 1 ,... ! ^),prW($ 7 , 2 (C/))(S 7?2 (^,...,^))). 
this implies pr^<I>(T/i, pr^$(?y2, z)) = pr( n )<I>(?7i + 772 , 5) . 



□ 



As in the classical case we therefore have (using the notations from Proposition 
EH): 

4.14 Proposition Let <3? be a projectable generalized group action on W x R q 
such that pr( n )<3? is a symmetry group of the algebraic equation A{z) = 0. Then <3? 
is a symmetry group of p3j. 



Proof. If U G G 7 

Thus 



is a solution of (34) then T-p T ( n )jj C 5a by Proposition 4.11 



■ pr(™)(#^t/) 



pr^%(T pT(n)u ) C 5 A , 



so that, again from Proposition 4.11, the claim follows. 



□ 



4.15 Definition Let X be a Q -complete generalized vector field. The n-th pro- 
longation of X is defined as the infinitesimal generator of the n-th prolongation of 
the generalized group action $ corresponding to X: 



pr(»}jr|, = ± 
an 



provided that pr( n )$ is Q- complete as well. In this case, both X and are called 
Q-n-complete. 



From Theorem 4.7 and Proposition 4.14 we immediately conclude 



4.16 Theorem Under the assumptions of Proposition [(.1 4 , let $ be a Q-n- 
complete generalized group action on R p x R 9 with infinitesimal generator X such 
that the conditions of Theorem are satisfied for A andpi^Q. If 

pi^X(A)(z) =0 VI G (K p x TZ q ) {n) with A(z) = 0, 

then <J> is a symmetry group of (fg^J. □ 

In order to be able to apply the same algorithm as in classical Lie theory for the 
determination of the symmetry group of a generalized PDE, the final step is to 
verify that the formulas for prolongation of vector fields carry over to generalized 
vector fields. 
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4.17 Theorem Let 



X = (x, u) -> + X] 



j=l c*=l 



6e a Q-n-complete generalized vector field with corresponding projectable group ac- 
tion $ on (MP x R q ). Then 



a=l J 

where J = (ji, j^), 1 < jk < P /or 1 < /c < n and 

i=l i=l 

Proof. Using the machinery developed so far, this is an easy modification of the 



proof of the classical result (see [29|, Theorem 2.36). □ 

We may summarize the results of this section as follows: In order to determine the 
symmetries of a differential equation involving generalized functions, the algorithm 
(as in the classical case) is to make an ansatz for the infinitesimal generators, cal- 



culate their prolongations according to Theorem 4.17 and then use Theorem 4.16 
to determine the defining equations for the coefficient functions of the infinitesi- 
mal generators. The defining equations now yield PDEs in Q T . Any solution of 
these equations that defines a ^-n-complete generator will upon integration yield 
a symmetry group in Q T . 

4.18 Example Scalar conservation laws of the form 

u t + F{u)u x = (37) 

arise in the kinetic theory of traffic flow. Here u denotes the density, and the 
propagation velocity F may be a strictly decreasing function of u with one or 
more jumps. A typical case is a unimodal flux function (whose derivative is F) 
with a kink at its maximum, as supported by experimental data |p.5| . Convolution 
with a nonnegative mollifier (p £ ) ee j allows to interpret F as an element of Q T (M) 
which is invertible. Thus our theory of symmetry transformations for equations 
with generalized nonlinearities applies. The determining equations are 

(ft + Fip x = 
+ Fr t + rF t + tpF u - F£, + F 2 t x + £F X = 

with infinitesimal generator v = £(x, t)d x +r(x, t)dt + <p(x, t, u)d u . As a particular 
solution we obtain v = xtd x + t 2 dt + (F'(u))~ 1 (x — tF(u))d u . The corresponding 
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generalized group action can be calculated explicitly in Q T showing that if u is a 
G T -solution to (j37|) then so is 

(x, t) F~ l (r]x(l + r]ty l + F(u(x(l + Tft)' 1 , i(l + ^)^ 1 )(1 + r]^ 1 ) 

In particular, a constant state u is transformed into a generalized solution to (|37| ) 
which, depending on the shape of F, will generally be associated with a piecewise 
smooth function. 



4.19 Example The nonlinear d'Alembert-Hamilton system 

Utt "ra Uyy u Z z — F(u) / qq\ 

uj - u 2 x - u 2 y - u 2 z = G{u) 

arises in the study of relativistic field equations Q and as a constraint in reducing 
the nonlinear wave equation to an ODE |l2|, |l3| . One of its symmetries is generated 
by the vector field v = ip(u)d u where the function <p has to satisfy 

Fip u - tpF u + Gifuu = 
2G(p u - t°G u = . 

In particular, in the isotropic case F = G = the function ip is arbitrary. In 
our theory it may be taken in £/ T (R) subject to the ^-completeness conditions 
formulated above. As an example of the possible behavior of generalized trans- 
formations, consider the vectorfield v = ip(u)d u where (p € ^ T (M) is the class of 
(ip £ ) £( zi with (p £ (u) = tanh(|). Thus (f(u) is associated with the jump function 
— sgn(it). Starting with a classical smooth solution u = u(x,t) G Oc(R 4 ) of the 
isotropic d'Alembert-Hamilton system ((|38|) with F = G = 0), the generalized 
symmetry transform generated by the vector field v turns u(x, t) into the general- 
ized solution U G t/ T (R 4 ) with representative 



: (x,t) = e Arsinh ^e v/£ sinh u ^ x,t ^ \ 



When rj > 0, it is straightforward to check that U is associated with the piecewise 
smooth function v(x,t) = u(x,t) + rj sgn(u(x , t)) . The generalized symmetry this 
way transforms smooth solutions into discontinuous solutions. 
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